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I.  INTRODUCTION 


The  theory  of  queues  is  rich  with  information  on  the 
classical  type  of  queueing  problems.  In  most  cases,  the 
solution  to  these  problems  is  straightforward  and  agrees 
with  one's  intuitions  concerning  actual  physical  situations. 
Quite  often,  however,  solutions  rely  upon  the  fact  that  ar¬ 
rival  and  service  distributions  are  independent.  In  n-server 
queues,  one  usually  stipulates  that  arrivals  enter  each  queue 
independently . 

The  purpose  of  this  thesis  is  to  investigate  three  types 
of  service  systems  in  which  arriving  customers  generate  jobs 
to  two  service  systems  simultaneously.  Thus,  while  the  serv¬ 
ice  systems  work  independently  of  each  other,  they  both  receive 
customers  via  the  same  arrival  mechanism. 

The  first  system  to  be  studied  is  the  two  channel  queue, 
where  each  service  channel  has  an  infinite  number  of  servers. 
The  second  system  is  a  modification  of  the  first  whereby  each 
service  channel  has  only  one  server  and  queues  are  allowed  to 
build  up.  Finally,  the  third  system  further  modifies  the 
situation  by  assuming  that  the  entire  system  has  a  maximum 
capacity  of  K  customers. 

The  following  chapters  will  show  that  the  infinite  server 
system  is  easily  solved.  However,  as  soon  as  one  restricts 
the  number  of  servers  to  one  per  channel,  several  difficul¬ 
ties  arise.  At  this  time,  only  partial  information  has  been 
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established  concerning  the  probabalistic  nature  of  these 
two  later  systems.  Hopefully,  further  investigation  can 
be  performed  utilizing  the  current  findings,  limited  as  they 
are . 
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II.  THE  TWO  CHANNEL  IKFiniTE  SERVER  QUEUE 


In  this  first  case,  it  is  assiomed  that  customer  arrivals 
occur  in  accordance  with  a  Poisson  Process  of  intensity  X  . 
Each  arrival  requires  two  separate  services  which  are  performed 
simultaneously.  The  two  service  channels  operate  independently 
and  each  have  an  infinite  or  unlimited  number  of  servers.  The 
servers  in  channel  one  work  at  an  exponential  rate  with  mean 


while  those  in  channel  t’^/o  work  at  an  exponential  rate 


with  mean 


Let  N  be  the  number  of  customers  in  the  system  at  equilib¬ 
rium.  Let  and  'be  the  number  of  jobs  being  performed  in 
each  channel  at  equilibrium. 

The  diagram  below  shows  the  system  for  an  arbitrary 
instant  in  time . 


\ 


M, 

Channel  1 


AA.  ^ 

Channel  2 


9 


Note  that  N  =  max  and  that  N^  and  N^  are  depend¬ 

ent  random  variables.  Furthermore,  because  the  service  chan¬ 
nels  are  independent  the  marginal  distributions  for  N^  and  N^ 
are  both  Poisson. 


In  order  to  describe  the  characteristics  of  this  queueing 
system,  it  will  be  necessary  to  apply  the  theory  of  the  filtered 
Poisson  Process.  Parzen  defines  a  filtered  Poisson  Process  as 
a  stochastic  process  which  can  be  represented  as 

A(t) 

X(t)  =  w{t,-Cn,Y„) 

where  A(t)  is  a  Poisson  Process  with  intensity  X  • 
is  a  sequence  of  independent  random  variables 
identically  distributed  as  Y  and  independent  of  A( t) . 

w(t,‘Xj^fY^)  is  a  response  function  and  is  the  time 

of  the  n^  arrival.  Here  w  indicates  whether  an 
arrival  is  still  present  in  the  system  at  time  t. 

In  the  context  of  queueing  theory,  Y^  will  be  the  sequence 
of  service  time  random  variables.  The  response  function  will 
be  defined  as 
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1  iff  t-t 

^  ’  '“n’  n-'  0  iff  Y  t-t 

n  n 


In  effect,  w  "counts"  the  number  of  customers  present 
at  time  t.  Thus  the  number  of  customers  present  can  be 
viewed  as 


A(  t) 


N(t)  = 


n=0 


w(t,Xn.  V 


A.  STEADY  STATE  DISTRIBUTION  FOR  THE  NUMBER  OF  CUSTOMERS 
IN  THE  SYSTEM 


Since  both  services  commence  simulaneously ,  the  service 
time  for  each  customer  is  Y^  =  max^S^.S^"^  where  S^'^Expo- 
nential  and  S^ Exponential  (  ^2 

ent.  Thus  for  all  n, 


P(Y^<t)  =  P(S^<t)  P(S2<t) 

-U,  t  -iut 

P(Y^<t)  =  (1  -  e  ^  )(1  -  e  ^  ) 


Parzen  shows  that  under  the  above  conditions  N(t)  is  Poisson 
distributed  with  mean  X  ^  (l-Fy(s))  ds.  Thus,  N(t)  has  mean 


-9, s  -9^s  -(On  +  0o)s 


+  e 


-  e 


^  ^  )ds.  In  the  limit,  N  is 


Poisson  distributed  with  mean  E(Y)  where  E(Y)  - 


Ml  +^2 


After  some  algebraic  simplification,  the  mean  is 
2  2  \ 

XU  +  Ml  +  M2) 

E(N)  =  - ^ - - - 

Ml  ^2  +M2) 


.  One  notes  that  E(N)  is 
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less  than  the  sum  of  E(N^)  and  ECN^) ,  which  is  to  be 
expected. 


B.  JOINT  MOMENT  GENERATING  FUNCTION  FOR  N^  AND  N^ 


Given  n  arrivals  in  (0,t),  the  arrival  times  are  distri¬ 
buted  as  n  order  statistics  from  a  uniform  (0,t)  random  var¬ 
iable  . 

Since  the  order  of  arrival  has  no  impact  on  service  (infi¬ 
nite  server  queue),  the  probability  that  an  arriving  customer 
will  still  be  present  at  time  t  is  merely  a  fimction  of  the 
service  distribution  and  the  time  of  arrival. 

Consider  an  arbitrary  arrival,  j,  in  (0,t).  The  arrival 
generates  one  job  for  each  service  channel.  Let 
YzjCt)  be  a  two  dimensional,  bivariate  Bernoulli  random 
variable  where  the  event  =  1  is  defined  as  the  j  job 

in  channel  1  being  present  at  time  t.  ~  ^  defined 

as  the  event  where  the  j  job  is  present  in  channel  2. 

Now  for  the  arbitrary  arrival,  let  x  be  arrival  time. 

Thus  X'v-' Uniform  (0,t)  and  the  joint  moment  generating 
function  for 


0  0 

(9^,02)  =  (p-]_  +  *^2® 


where  p^  =  1  -  F^(t-x)  ,  p^  +  q^  =  1 


p^  =  1  -  F2(t-x)  ,  P2  +  q2  =  1 


12 


and  and  F^  are  the  service  time  distribution  functions  for 
channels  1  and  2  respectively. 

Thus,  unconditioning  on  the  arrival  time  X, 


Y,.(t).Y2.(t) 


^®1’®2^  =  ^  (l-F^(t-x)  +  F^(t-x)e^^). 


(1  -  F„(t-x)  +  F_(t-y)e  ^)-dx 


Letting  w=t-x  for  a  change  of  integration  variable 


Yi^.(  t)  »Y2^.(t) 


0. 


(w)e  -")■ 


(l-F^Cw)  +  F2(w)e  dw 


Considering  all  n  arrivals  in  (0,t),  each  arrival 
creates  an  independent  bivariate  Bernoulli  random  variable. 
Hence,  letting 


n 


n 


X  (t)  =  ^  .  . 

j=l  ^ 


Y-,  ,(t) 


X3(t)  =  _ 


Yo.(t) 


then  the  joint  moment  generating  function  for  X^(t)  and  X2( t) , 
conditional  on  n  arrivals  in  (0,t),  is 


X^(t).  X^Ct) 


{e^.ep  = 


[■ 


Ylj(t),  Y2j(t) 


(6^, 83) 


n 


r- 1 


(0^,02)  =  ^  ^  (l-F^(w)  +  F^(w)e 


0. 


X^(t),  X^Ct) 


(1-F,(w)  +  F„(w)e  ‘^)dw 


n 
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Thus  letting  N^(t),  N2(‘t)  the  numbers  present  at  time 
t  for  any  number  of  arrivals  in  (0,t) 


IVI 

N^(t),  N2(t) 


( 61.62^ 


n=0 


X^(t) 


X2(t) 


(©1,62^ 


It  follows  that, 


M  (©T  ,0.)  =  e"  e^ 

N^(t),  N^Ct)  ^  ^ 


where  A 


0 


0 


( 1-Fj^(  w)+F^(  w)  e 


■)  ( l-F^C  w)+F2(  w)  e  '^)dw 


Letting  and  be  Bernoulli  random  variables,  specific- 
cally 

Y^'^  Bernoulli  (p^^)  where  p^^  =  l-F^(w) 

Y^'^  Bernoulli  (P2)  where  p^  =  l-F^Cw) , 

then  the  moment  generating  function  can  be  expressed  as 


M 

Ni(t)  ,  ^^it) 


vl  0-,  +  0  0-|  6p 

A|_a(e  ^  ^  -1)  +  b(e  -l)  +  c(  e  -1) 
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where 


a  = 


P1P2  ‘iw 


t 

0  =  ^  dw 

0 


Now  consider  the  limiting  form  of  the  above  function  as  t 
First 


lim  a  = 
t-^”« 


-H^w  -/UpW 
e  e  dw 


1 


lim  b  = 

t— ^ 


-^w  -  U^w 

e  (1-e  )  dw 


_1 _  _  1 

^1  ^1  '*■'^2 


M2 


lim  c  = 


-M  W  -Mn  W 

e  (1-e  )  dw 


1 _ 

^2 


1  ^  ^1 _ 

+x(2 


Then 


M  (9,  ,0.)  =  lim  M  ^®1’®2^ 

N^.N^  t-^'O  N^(  t)  .N^Ct) 


X 


e 


XA  T  +  -U. 


( e 


V®2 


-  1)  + 


M 


XI  ^^Tu^Tirp’ 


-  1) 


M  1 _ (^2  _ 

^2(^1  +^2) 
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It  should  be  noted  that  the  marginal  generating  functions, 
and  hence,  N^,  N2  are  Poisson. 


\  0-, 


=  e 


/C^^(e  1) 


Now  consider  the  second  mixed  moment,  E  [ninJ  . 


E  =  li™ 


A 


Sj;^  0  ^®2 

92-»0 


in  


=  \ 


-1  *-”2 


®r®2 

e  ^  + 


and 


M 

)-  Q. 


1 

-L  Q  ^  + 


-Ul+  /U2 


'^2  ®l'\ 

j  ^92 


X  er«2 


"1  *^z 


'*'Ni,N2^®1'®2* 


Therefore 

E 


[niN^]  = 


X 


lim 


AK 


1  0 


X 


Since  N^  Poisson 

1  M 


M 


\ 


^8; 


"1  ■  ^2 


and  N„  Poisson 


X 


16 


then  lim 
®2-"0 


M 


^9. 


=  E  U.l  = 


X 


and 


E  CNj]  = 


\ 


M 


Thus,  the  correlation  between  and  is 


E(N^N2)  -  E(N^)  ecn^) 


V(N^)  ^(N^) 


J- 


1 


v^2 


**1^2 


X 


X  ,  X 


/U  '  ■«. 


(  -7^) 


X  X 


Ax. 


1 


J 


1  'Cl  p 

-7^ 


,  which  is  independent  of 
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III.  THE  TWO  CHANNEL  QUEUE  WITH  A 
SINGIE  SERVER  IN  EACH  QUEUE 


As  a  modification  to  the  first  system,  each  service  channel 
will  now  he  limited  to  a  single  server.  The  servers  in 
channels  one  and  two  work  with  exponential  rates. and 

respectively.  As  before,  N  is  the  number  of  customers 
in  the  system  at  steady  state.  and  are  the  nimnber  of 

jobs  to  be  done  in  each  channel.  The  jobs  are  served  using 
a  first  in- first  out  discipline.  The  arrivals  are  generated 
by  a  Poisson  Process  of  intenstiy  X  • 

A.  STEADY  STATE  BALANCE  EQUATIONS  FOR  THE  SYSTEM 


Since  the  system  has  both  exponential  interarrival  and 
exponential  service  times,  the  balance  equations  for  the  two- 
dimensional  steady  state  probabilities  can  be  expressed  as 


follows.  Letting  P  =  P(N-=n, ,  N„=n«)  one  has  that 

>  n  2  -L  -L  ^  ^ 


1<  n^ 

<  ‘O 


(  X+M^) 


n^,  0 


=  1.0 


n^,l 


1  4n^  < 


(X 


"2> 


0  ,n. 


=  1X 


P 


0 ,  n 


+  1  + 
2 


V 

1 ,  n2 
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1  ^ 


X 


00 


=  Xl 


1  10 


+  p 


2  01 


The  transition's  between  states  can  best  be  visualized 
using  the  following  transition  diagram 


One  might  suppose  that  the  balance  equation  could  be  solved 
using  difference  equation  techniques .  This  does  not  seem  to 
be  the  case.  First,  most  techniques  used  to  solve  partial 
difference  equations  of  the  above  type  assume  a  solution  of 

the  form, 
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Under  the  above  assumption,  one  could  solve  for  P 
and  express  the  result  as 


K 

P  =  ^ 

n,  ,n,  Z— 

^  ^  j=l 


^1  *^2 

c.  a.  b.  where  (a.,b.) 
J  J  J  J  J 


j  =  1,  2,  .  .  .  K 

are  the  K  roots  of  the  characteristic  equation  and  c^'s  are 


arbitrary  constants.  However,  since  P 


is  a  function  of 


two  dependent  random  variables  and  N^,  it  is  known  that 


'^1*^2 


cannot  be  factored  into  a  separable  product  similar 


ni  n 

to  a  b  .  If  this  were  the  case,  then  and  would  be 
independent . 

Therefore,  for  each  probability  there  are  Vmo  or  more 

arbitrary  constants  to  be  evaluated.  If  one  solves  the  major 

difference  equation  for  the  interior  points  of  the  first 

ni  n^ 

quadrant,  i.e.  assume  P^  ^  =  a  b 


+  ^^n^-l,n2-l  l^n^+l,n2  ^  ^2^n3_,n2+l 


the  roots  are  a=l  b=l,  a=l  f=- 


M., 


and  a  = 


\ 


A 


b=l.  Thus  the  general  solution  is 


^  ^1  n2 

’n,,n,  =  2-  <=j(ni.n2) 

^  j=l 
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or 


rii.n^ 


'X  ^2  \  n 

CiCrii.n^)  +  ^  '*’ 

^  -  1 


Even  using  all  of  the  boundary  conditions,  one  cannot 

find  a  particular  solution  to  the  above  equation. 

One  point  to  note  is  that  P  can  be  expressed  as 

■^1^2 


^1^2 


=  P  +  P  -  P(n,  v;  n^) 


n 


n 


Thus  if  one  assumes  that  ^2^^1’^2^ 
actually  (l  -  )  and  (l  -  -jj- — )  respectively,  the 


'2 

particular  solution  takes  the  following  form 

\,n3  =  ^ 


n- 


) 


Again,  using  the  boundary  conditions  only  results  in  a 
system  of  equations  for  C(n^,n2)  which  is  equivalent  to  the 
original  set  for  P 

It  appears  as  though  the  boundary  conditions  are  of  a 
non-separable  type,  and  consequently  there  is  no  explicit 
solution  for  the  balance  equations. 

B.  SOLUTION  TO  THE  BALANCE  EQUATIONS  USING  THE  PROBABILITY 
GENERATING  FUNCTION 

The  joint  probability  generating  function  can  be  expressed 
as 

G(s,t)  = 

ni=0 
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Using  the  balance  equations  one  can  derive  the  following 
expression  for  the  joint  generating  function. 


G(  s , t) 


(  \+  ^  2^  ^  s^t^ 


s«2 


In  other  words,  G(s,t)  is  expressed  as  a  function  of 
s,  t,  and  the  boundary  probabilities.  If  one  knew  these 
probabilities,  G(s,t)  could  be  solved  in  closed  form.  However 
such  is  not  the  case.  It  should  be  pointed  out  that  the  same 
type  of  equational  fonn  occurs  when  finding  a  generating 
function  for  queues  which  have  either  Gamma  interarrival  or 
Gamma  service  times.  However,  the  boundary  probabilities 
are  finite  in  number  and  can  be  solved  for  using  Rouche's 
theorems  on  the  roots  of  complex  functions.  It  appears  as 
though  this  method  is  fruitless  for  the  above  function. 

As  in  the  infinite  server  queue  of  Chapter  2,  the  marginal 
distributions,  ^  (n^)  and  known.  They  are 

geometric  and  can  be  expressed  as 


n^  =  0,  1,  2 


0,  1,  2 


n 


2 


As  one  would  expect,  the  marginal  generating  functions 
are  also  geometric. 
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(t-l)Xl  ^n,,0 


G(l,t)  = 


n^=0 


G(l,t)  = 


(  t-  Xt^-t/U^-  Ug 

(t-1)  (1-  )  (t-l)M  2  (1- 

2  2 


-\t^+(  X  +>U2)-t-Xl2 


-  Xt(t-l) 


"2(1  -  A 


X 


1  - 


X 


JU 


G(l,t)  = 


-  X  "fc  +  ^ , 


1  - 


IX 

Mr. 


which  is  the  generating  function  for  a  geometric  random 
variable  with  parameter  1  - 


1  - 


A 

M  1 


Likewise  G(s,l)  = 


Thus,  little  has  been  gained  using  the  generating  function 
approach. 


C.  AN  ATTEMPT  TO  FIND  THE  JOINT  GENERATING  FUNCTION  BY 
CONDITIONING  ON  ARRIVING  CUSTOMERS 


In  Chapter  2,  the  joint  generating  function  of  N^  and  N^ 

was  found  using  a  filtered  Poisson  Process.  In  essence,  the 
method  worked  because  the  stochastic  dependence  of  N^(t)  and 
N^Ct)  was  handled  by  conditioning  on  both  the  number  and 
arrival  times  of  n  customers,  given  an  arbitrary  time  t  since 
the  process  started.  Thus  N^(t)  and  N^Ct),  conditioned  on 
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the  above  events,  only  depended  upon  their  respective  service 
distributions . 

A  possible  approach  to  finding  the  joint  moment  generating 
function  for  N^(t)  and  N2(t)  is  to  handle  the  stochastic 
dependence  by  the  same  conditioning.  Thus,  given  the  number 
of  arrivals  at  time  t. 


n=0 


Letting 


P  [A(t)  =  n] 
^  =  E  [e  1  ^  /A(t)  =  nj 


it  follows  that 


t  t  t 


0  = 


E  [e  ^  /A(t)  -  n,  S^,S2,  .  .  .S 


°  ^1  Vi 


n 


—  dsT  dSo  • • • ds^ 
n  1  2  n 


r  }  \  \  1  ' 

®  I®  /N(t)  =  "iSiSi-’-Sn  J 


S,  S  T 
1  n-1 


'H2N  (t) 

_e  ^  /N(t)  =  n^S^, 


.S 


n  ^ 


J  ds.  ...  ds 


,n  "*"1 


n 
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Now  the  problem  is  to  find  the  moment  generating  function 
for  N^(t)  and  N^Ct)  conditioned  on  the  two  given  events.  It 
appears  that  this  is  not  a  trivial  problem. 

One  might  recall  that  in  the  infinite  server  queues,  N^(t) 
and  N^Ct),  conditioned  on  the  same  events,  were  actually 
binomially  distributed  with 

p  =  i  ^ -  dt 

0  t 


The  probability  that  any  customer  would  be  present  at  t  was 
independent  of  the  presence  of  other  customers.  However, 
with  a  waiting  line  in  each  queue,  a  dependence  is  generated. 

If  the  arrival  is  present  at  time  t,  so  are  the  arrivals 

j+1,  j+2,  ...n. 

To  this  date,  the  conditional  distribution  of  N^(t)  (or 
N2(t))has  not  been  determined.  One  point  to  consider  is  that 
whatever  the  conditional  distribution  is,  it  must  be  such 
that  the  limiting  distribution  of  N^( t)  (or  N2(t))  is  geometric. 

It  should  also  be  observed  that  the  established  theorems 
concerning  a  filtered  Poisson  Process  do  not  apply  when 
waiting  lines  develop.  Even  though  N^(t)  can  be  expressed  as 


the  response  functions  are  correlated.  If  w^ 

4“  Vi 

customer  then 

w^  =  1  for  ^  -c  i  ^  n. 


1  for  the 
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D.  ANALYSIS  OP  ABSORPTION  PROBABILITIES  USING  RANDOM 
WALK  THEORY 


Consider  the  situation  where  the  state  of  the  system  is 
given.  That  is,  at  the  present  time,  N^  =  i  =  j.  What 

is  the  probability  that  the  system  will  reach  one  section  of 
the  boundary  before  the  other  section? 

The  attempt  to  answer  this  question  used  the  theory  of 
a  random  walk.  The  problem  reduces  to  solving  an  infinite 
system  of  difference  equations  of  the  form 


P 

n 


1 


.n 


2 


\ 

X  +  U, 


u 


n^-l. 


n2-l 


'  n^+ 1 , 


n< 


*  X  +  Ui+^2 

for  l^n^-^®®  l‘^:n2'^‘=^ 

In  this  case,  the  boundary  is  partitioned  into  two  classes, 
following  an  approach  introduced  by  Feller  in  An  Introduction 
to  Probability  Theory  and  Its  Applications,  Volume  I.  Using 
this  method  of  analysis,  the  boundary  points  are  either  given 
values  of  1  or  0.  Thus  if  P  o  probability 

f  V  u ,  112 

of  being  absorbed  by  this  part  of  the  boundary  is  1.  Other¬ 
wise,  the  absorption  will  take  place  at  the  other  boundary 
points . 

Unfortunately,  the  solution  to  this  system  is  as  obscure 
as  the  original  system  of  balance  equations  explored  earlier. 
At  this  point,  no  solution  has  been  found. 
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E.  SOME  ADDITIONAL  FACTS  CONCERNING  THE  TRANSIENT  DISTRIBUTION 
OF  WAITING  TIMES  FOR  ARRIVING  CUSTOMERS 


The  following  mterial  uses  the  queueing  theory  approach 
espoused  hy  Feller  in  An  Introduction  to  Probability  Theory 
and  Its  Applications.  Volume  II.  The  intent  is  to  develop  a 
closed  form  expression  for  the  waiting  time  of  the  j  arrival, 
given  that  n  customers  had  arrived  at  time  t. 

It  will  he  noted  that  the  theoretical  results  are  highly 
intractable.  Perhaps  further  research  will  simplify  the 
following  relationships. 

Let, 


X  =  B  -  A 
n  n-1  n 


=  •••  ""n 


n-1 


n 


n 


A. 


* 

=  B  -  A 


=  Z  B,- 

i=0  j=l 

n  =  number  of  customers  -  (start  with  customer  O) 

A  =  time  between  arrival  of  n-1^^  and  n^^  customer 
n 


B  = 
n 

service  time 

of  n*>' 

customer. 

W  = 
n 

waiting  time 

of  n*'' 

customer 

W  = 
n 

max  VO,  S^, 

S^,  . 

••• 

If  S^  is.  max  ^0,  S^, 

®2’  ■ 

•••  ^n] 

then  S^>  0, 

S^>  s^,  > 

Sj.  .  . 

• •  ’  ^k-l’ 

S  < 
k+1 

^k’  ^k’ 

•  •  •  • 

Bn-1  ®n 
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The  event  0,  Sj^  >  . Sj^  >  Sj^_. 

independent  of  4  Sj^,  S^2<  S,^ . <  S^? 


Also  the  event 


I 


^  ^  ^1’  •  •  •  •  Sj^  > 


is  equivalent  to 


Si  >  0,  >  0 .  ^k  ° 


the  event 


l5+l 


"  Sv  <  0,  Sj^2  -  <  0,  ...  S^_T  <  0,  -  S,,<.  0 


n-1 


n 


is  equivalent  to 

^  <  0,  <  0,  S^<  0 . S 


n-k-1 


<  0, 


Sn-k<  0 


Thus,  is  max  ^0,  S^, 


S„ . S 

2  n 


if  the  two  independent  events 


Si>0,  S2>0 . Sj^>0?  S2<0 . ^n-k"^  ^ 


are  true . 


Let  Pk“^/^i'^°'  ^2^°’  ■■*'  ^k^  .  .  .  n 


‘In-k  PjS^^O,  Sg  <  0 . ^n-k^°(  ^0  ^  ^0  ^ 


Then  P(Sj^  is  max)  =  q 


k  ^n-k 


Thus  the  distribution  of  =  max  -^0,  S^,  ....  S^ 

n 


is  P(W  <t)  =  ^  p(s.  <  t)  P(S.  is  max:) 
i=l  ^  ^ 


or  W  =  0  if  0  is  max. 
n 
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Therefore 


n 


P(W  <  t)  =  ^  P(S.  t)  p.  q  . 

n  Z. 1  -^1  ^n-i 

i  =1  - 


P(W  =  0)  =  q 
n 


k 

Consider  (B^_^  -  A^)  = 

i-1 


p(s^>o)  =  p(b;.^>a;  ) 

oo 

Gamma  (K^, 

■»• 

Aj^  ^  Gamma  (K^, 

>U;l) 

X) 

p(s^>  0)  =  t) 

a 

-f  *(t)  dt 

h. 

^  k-1 

=  \  Z 

J  j=o 

(X4,t)  X 

j!  (k-1) ! 

CD 

1 

>" 

c+ 

dt 

U  T 

=Z 

j=0 

x'^ 

( ic^-  i‘h !  . 

(k-1) ! 

(  X+A^)‘"+j 

P(Sk>0) 
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Recognizing  this  to  be  the  CDF  of  a  type  of  Negative 


Binomial  distribution, 

=  P(2^<K)  where  Z^Negative  Binomial  with 

*  ,  v  . 

parameters  P  =  ( — r - ) 

A  +U^ 

and  K 


3 


Z  =  0,  1,  2, 


Note  also  that  2^  can  be  expressed  as  the  sum  of  independent 
geometric  random  variables. 

X.  (X)  = 

i  X  +M.  V  X  + 


2,  =  X^+X  + 
k  12 


•  •  •  • 


where 


4 


X  =  0,  1,  2,  _ 

By  symmetry,  the  event  P(S^<0) 
can  be  expressed  as  )  “  ^^^k  ^  ®k-l^ 

thus  P(Sj^<0)  =  PCW^^K) 

*  ^1 

where  Negative  Binomial  with  parameters  P  =/ -  and  K. 

^  A  +M^ 

can  be  expressed  as  the  sum  of  K  independent  geometric 

random  variables. 
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=  Y^+Y2+  .  .  . 


where 


Y  (y)  = 


\ 


y 


Y. 

1 


X  +  ^ . 


X  +  . 


Now  given  a  sequence 


y  =  0,  1,  2, 


^ ’  ^1^2'  • * ' * 

consider  computing  is  max) 

P(S^  is  max)  = 

=  P(S^>  0,  0,  _ ^  ^  P(z^<  1,  <C2 . ^k 


Pk  "  X^+X2<2,  X^+X2+X^<3, 


X1+X2+X3+.  .  .  .X^<K) 


The  last  event  can  he  described  as  a  sequential  modified  random 
walk  where  X^  =  0,  1,  2,  ....  and  the  event  can  only  occur  when 
each  partial  sum  is  less  than  K. 

As  an  example: 

X^  must  =  0 


X^  can  equal  0  or  1 
since  X^  +  X^  <  2 


etc. 


31 


A  pictorial  representation  of  the  above  sequence  is  shown 


below  for  K  =  5 


This  sequence  satifies  the  event 


n 


<  4 
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In  order  to  calculate  the  probability  of  the  event  K=5,  i,e. 

=  P(X^^1,  X^+X2<2,  - X^+X^+X^+X^+X^  <  5) 

one  must  calculate  the  probability  for  each  path  starting  at 
(1,0)  and  terminating  at  either  (5»0)  (5»l)  (5,2)  (5,3)  while 
remaining  in  the  acceptable  region  of  transitions. 

It  so  happens  that  there  is  a  pattern  to  these  paths. 

In  order  to  terminate  at  (5»0) 
all  X^  =  0  P((1,0)-^  (5,0))  =c<^ 

where°<.=  r - 

X  + 

In  order  to  terminate  at  (5»1)  all  X^  =  0  except  for  one: 

P((l,0)— ^  (5,1))  =^^(1-®^)  One  must  multiply  this  probabil¬ 

ity  by  a  factor  reflecting  the  number  of  paths  from  ( 1 , O)— > (5 , 1) • 
This  will  be  considered  shortly. 

The  other  termination  probabilities  are 

(1,0)-^  (5.2)  :  oC^(l--C)^ 

(1,0)  (5»3)  !  ^<^(1-^^)^ 

(1,0)  (5,^)  J  ^<^(1- 

Thus  we  have  shown  that 

P^  =  P(X^  1,  X^+X2<2 . X^+X^+X^+X^+X^<  5) 

i=l 

K  ^ 

In  general  P^  =  i-1^  ( 1- ‘^) 

i=l 

Now  what  about  the  coefficients 
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One  can  analyze  the  possible  paths  to  a  given  termination 
point  by  an  induction  method.  Consider  the  cases  K=l,  2, 
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For  K=1  P^=P(X^^1)  =  P(X^=0)  =  -C 

For  K=2  the  two  termination  points  (2, 0)^(2, 1)  have  only 
one  path  to  each  of  them  hence  C^q  =  1  ^21  ~  ^ 

P^  +  X^(l-0 

For  K=3,  the  three  termination  points  are  (3.0)  (3,1)  (3,2). 

Consider  (3.0).  It  can  only  he  reached  from  (2,0). 

Since  the  number  of  ways  from  (1,0)-^  (2, 0)  =  1,  then 

^30  ^  ^20 

Consider  (3.1) <  It  can  be  reached  from  (2,0)  and  2,1). 
Hence  C^^  =  C^q  +  ^21'  likewise  C^^  “  ^20  ^21’ 

Thus  it  appears  that  one  can  compute  the  using  a 

modified  form  of  Pascal's  triangle. 

K  n=7 

1  1 

2  11 


3 

1 

2 

2 

k 

1 

3 

5 

5 

5 

1 

k 

9 

14 

14 

6 

1 

5 

14 

28 

42 

42 

7 

1 

6 

20 

48 

90 

132  132 

etc . 

where  each  row  starts  with  1  and  each  successive  member  is 
equal  to  the  sum  of  all  row  members  in  the  preceding  row 
starting  from  the  left  up  to  but  not  including  the  same 

column. 
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i.  e. 


"K 


i  =  %  °K-1,  j 


Now 


%-k  ^2^  ° . ^n-k  ^ 

=  P(Y^<  1,  Y2<  2,  .... 


*  1 

where  the  Y. Geometric  (P  =  r - - —  ) 

^  X  + 


By  analogy 


n-k 

^n-k  ”  ^n-k,j-l  P 

j=l 


n-k 


Hence  P(S,  is  max)  =  p,  q 


k^n-k 


il-f) 


j-1 


? 


=  ( 


jA 


X  +M. 


n-k 


P(S^  is  max)  =2^^  ^k,i-l^  (1-0 


i-1 


j=l 


'n-k, j-1 


*  "■’'(i-p ) 


3-1 


To  date,  no  method  has  been  foimd  to  simplify  the  above 
expression  because  of  the  complexity  of  the  coefficients,  c. . 

o 
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IV.  THE  TWO  CHAMEL  SERVER  QUEUE  WITH  SINGLE  SERVERS 

AND  LIMITED  QUEUE  CAPACITY 


In  this  chapter,  the  system  will  be  further  modified  by 
placing  a  limit  on  the  total  number  of  customers  allowed  in 
the  system.  Thus  if  N=K,  this  implies  that  N^  and  N^,  the 
number  of  jobs  in  each  queue,  are  also  limited  by  K.  Since 
N  =  max  then  the  ordered  pair  (N^,N2)  are  restricted 

to  states  bounded  by  a  square  region  in  the  first  quadrant. 
The  diagram  below  illustrates  the  aLLlowed  transitions  for  K 
equal  to  3. 
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The  balance  equations  for  the  K  capacity  system  are 


complicated,  but  finite  in  number. 


(X 


+  UU+  =  X  P  ,  ,  +  n.p  +  K  P 

1.^  2^  '^l'^2  n^-l,n2-l  1  n^+l,n2  .  2  n^, 


1  K  1  £  <  K 


I4n^<  K 


(  X+ Ai2)Po^n2  “  '^2^0,n2+l  ■^^l^l,n. 


l-n2<K 

{XX  2)PK,n  "  ^  ^K-l,n2-l  “^^2  ^K,n2+1 


If  n2<K 


(ii,+  tt.)p  „  =  X  P^  ,  ^  ,  +X(  p 

'  1  2'  n^^tK  n^-l,A-l  i  n^^+l, 


K 


^  ^00  Ao  ^  '^2  ^01 


1  f-  K 


^1^K,0  '^2^K,1 


'^2^0,K 


K 


"  ^^K-1,K-1 
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There  has  been  a  great  deal  of  research  performed  in  the 
area  of  finding  solutions  to  large,  finite  state,  queueing 
systems.  The  purpose  of  this  work  was  to  avoid  the  use  of 
computer  simulation  or  standard  iteration  procedures.  Instead, 
an  investigation  into  the  special  structure  of  this  finite 
system  was  conducted.  Although  a  solution  was  not  obtained, 
several  interesting  facts  were  documented  for  use  in  further 
research. 

It  should  be  noted  that  the  K  capacity  system  is  more  diffi¬ 
cult  to  attack  because  the  marginal  probability  functions  for 
and  N2  are  not  known.  In  the  previous  chapters,  arrival 
dependency  did  not  affect  the  marginal  distributions.  Here, 
one  sees  that  when  either  service  channel  has  K  jobs,  arrivals 
are  turned  away. 

A.  ANALYSIS  OF  THE  TRANSITION  DIAGRAM  FROM  THE  VIEWPOINT  OF 

NETWORK  FLOW 

One  method  of  observing  patterns  among  the  steady  state 
probabilities  is  to  view  the  transition  diagram  as  a  network. 

By  slicing  the  network  with  a  cut,  one  can  equate  probability 
flow  across  the  cut.  This  is  actually  doing  no  more  than 
adding  together  several  balance  equations.  However,  the  method 
does  facilitate  a  rapid  development  of  relationships  between 
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the  probabilities.  A  transition  diagram  for  K  equal  to  two 
is  shown  below. 


N 


2 


By  equating  flow  across  the  cut  one  finds  that 

=«  1^22  ^V20 


This  equation  could  have  been  found  by  adding  together  the 
three  balance  equations  for  ^20’  ^21'  ^22  follows. 

^1^20  =  ^2^21 

'W^21=  ^\0^^2^22 

-h  ( L*''“2^  ^22  "  ^  ^11 

*"1^21  ^^1^22  =  *  ^^10 
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Now  let  us  consider  an  arbitrary  value  K.  For  any  Finite 
K,  the  same  type  of  cuts  can  be  performed.  If  one  writes  all 
the  equations  resulting  from  vertical  cuts  such  as  the  one 
above,  it  follows  that, 


K-1 

X  ^  p 


K 


n^=0  ^''^2  ^  1^0 


K-1 

X  X  ^ 


K 


l,n.,  '^1  ^  - 


n2=0 


2,  n. 


n2=0 


X 


K-1  K 


X  K-l,np  1 

n2=0  n2=0 


■K,n. 


Using  the  fact  that  P(N^  =  i)  = 


K 


n2=0 


^,n2 


0  4i^K,  one  sees  that  the  vertical  cuts  actually  show  relation¬ 
ships  between  the  marginal  probabilities.  In  particular, 

X  P(N^  =  0)  -XPq.K  "  ^1  ^(^1=1) 

Xp(N^  =  1)  -Xp^^j^  =  P(N^=2) 

X  P(N^=K-1)-XPj^_^^j^  P(N^=K) 


In  order  to  eliminate  the  P . •  terms ,  one  can  add  all  of 
the  above  equations  to  obtain 

\  [l  -  P(N^  =  K)j  -X[p(N2  =  K)  -  =lt^[l-P(N^  =  O)'] 
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Note  that  the  only  probability  which  is  not  a  marginal  is 


R 


K,K* 


In  an  analagous  fashion,  all  the  horizontal  cuts  in  a  K- 
level  transition  diagram  yield  the  following  equation. 


Eliminating  P^r  ^  from  both  equations  yields  the  surprising 
result  that 


At  the  least,  the  above  relationship  could  be  used  as 
a  check  against  solutions  obtained  from  iteration  or  simulation. 
The  above  results  also  show  that  the  marginal  probabilities 
may  be  more  efficiently  solved  for  if  one  knows  the  boundary 
probabilities  in  the  original  transition  network. 

Another  interesting  set  of  relationships  can  be  found  by 
combining  the  first  two  equations  found  from  the  vertical  cuts. 
Recall  that 


X  P(N^=1)  ^'^l 


Additionally  using  the  balance  equation  for  node  it 

follows  that  ^*^2  ^0  K  ^^1^1  K*  equations  can  be 

combined  to  yield 
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P(N^=0)  =  +  P(N^=1)  P(N^  =2) 


2 

^1 

It  follows  by  symmetry  that 


\  Qj 

M2 


P(N2=0) 


(\  +  a^)  P(N2=1)  -A2P(N2=2) 


Thus  one  has  obtained  two  equations  strictly  in  terms  of 
the  msirginal  probabilites , 

If  one  continued  this  pattern,  it  would  not  be  possible  to 
arrive  at  an  entire  system  of  pure  marginal  equations.  However, 
the  above  equations  can  be  used  as  check  against  any  analytical 
procedure  which  might  seek  solutions  to  the  marginal  probabili¬ 
ties.  This  check  is  used  later  in  this  paper  to  show  that  a 
proposed  semi-Markov  model  for  marginal  probabilities  is 
incomplete . 

Several  other  cuts  were  attempted,  but  the  results  always 
appeared  to  re-establish  the  previously  stated  equations. 


B.  REDUCTION  OF  THE  NUMBER  OF  BALANCE  EQUATIONS  BY 
SUBSTITUTION 

Consider  again  the  system  for  K  equal  to  two.  It  will  be 

shown  that  the  nine  homogeneous  balance  equations  can  be 

reduced  to  three  by  a  method  of  substitution.  In  terms  of 

2 

an  arbitrary  K ,  this  translates  to  reduction  of  (K+1) 
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"balance  equat:ions  into  a  system  of  2K-3  balance  equations  in 
2K-3  unlcnowns.  The  method  is  as  follows. 


In  theory,  it  is  possi"ble  to  express  each  probability 

in  the  set  circled  in  terms  of  boundary  probabilities  ,  , 

^2  2  ^12*  ^  sense,  one  is  '* walking”  from  one  boundary 

to  another  on  the  opposite  side  of  the  network.  In  finding  all 
P^^^s  in  terTns  of  2  and  P^  2*  only  balance 

equations  '/i^ich  are  not  utilized  are  those  for  P^  Pq  q, 
and  Q*  As  an  example,  Pq  ^  can  be  expressed  in  terms  of 
2  ^2  2  using  the  balance  equation  for  P^  2»  "that  is 


^1+^2 


M 


■0,1 


\ 


■1,2 


p 

X  2,2 


Thus  if  all  the  P.  .’s  are  expressed  in  terms  of  P^ 

xj  ^'*^2 

where  P^^  ^  where  15n^-^K,  and  one  can  use 

the  remaining  boundary  equations  for  P^  l£n^<-K,  Pq^^ 

1  2 

16-n^^K  and  P«  «  to  arrive  at  a  homogeneous  system  of  three 

2  U  >  U 

equations  in  three  unknowns.  For  the  above  case. 
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M 


0,2  ^2 


U 


0,1  ^  1,2  X  2.2 


M 


JU 


0,0 


2.2  X  2.1  n.2 


1,0  -  ^ 
A  o 


iA 


2,1 


X  2,2 


2,0  ju  ^  2,1 


.+  A„ 

^2,2 


Now  if  one  writes  the  balance  equations  for  nodes  P 


0,1' 


Pq  q,  and  Q  using  the  above  equations,  the  result  will  be 
three  homogeneous  equations  in  three  unknowns,  P,,  ,  ,  P_ 

^  f  J.  <L  ^  ^ 


and 

It  must  be  pointed  out  that  there  did  not  appear  to  be  any 

obvious  pattern  to  the  build  up  of  coefficients  as  one  "walks” 

through  the  network.  All  attempts  made  to  find  closed  form 

expressions  for  the  coefficients  were  to  no  avail.  In  theory, 

each  P- .  can  be  expressed  as 
^  J 


K-1  ^^1 

’i'j  "  ^K,m  ^  ^n,K  ^  ^K,K 
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The  above  substitution  is  no  different  from  a  Gaussian 
elimination  scheme  that  could  be  used  for  any  system  of 
equations.  However,  it  was  hoped  that  these  equations 
contained  special  structure  and  would  lend  themselves  to  an 
analytical  solution. 

C.  COMMENTS  CONCERNING  THE  MATRIX  STRUCTURE  OF  THE  BALANCE 
EQUATIONS 

One  can  write  all  the  balance  equations  in  the  form  AX=B 

2  2  p 

where  A  is  (K+1)  by  (K+l)  ,X  is  the  (K+l)^byl  column  vector 

2 

for  the  ^ '  s  and  B  is  a  (K+l)  by  1  column  vector  which  has 
all  but  one  element  equal  to  0.  The  result  is  a  system  of 
(K+l)  -1  homogeneous  equations  and  1  normality  equation 


Z—  P,"  =  1.  The  normality  equation  may  replace  any  one  of 

ij  ^ 


2 

the  original  (K+l)  homogeneous  equations  and  thus  insure  that 


the  non-homogenous  system  has  a  unique  solution. 

After  a  great  deal  of  experimentation  with  reordering  the 
rows  and  columns  of  the  coefficient  matrix,  it  was  concluded 
that  there  existed  no  special  structure  in  the  equation  set. 
At  first  glance,  one  might  be  attracted  to  the  regularity  of 
the  coefficient  matrix.  As  an  example  writing  the  balance 

equations  for  K=2  in  descending  order  of  the  second  and  then 
first  coordinate  one  arrives  at  the  following  system. 
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I 


-X 

^22 

0 

— M-U 

"X 

H 

c\J 

0 

^20 

0 

-Uj 

' 

^12 

0 

-X 

hi 

- 

0 

^10 

0 

C\J 

o 

0 

Xv>Ux 

^01 

0 

) 

\  \  1 

1  '  .  \ 

»  1 

1 

o 

o 

1 

_ 

' 

i 

••  ^ 

The  coefficient  matrix  is  strongly  diagonal  in  structure. 
However,  no  special  method  was  found  to  solve  the  above  system 
without  resorting  to  the  standard  elimination  schemes  already 
known..  In  all  cases,  it  was  necessary  to  completely  triangular- 
ize  the  matrix  and  then  use  back  substitution  to  solve  for  the 


Originally  the  intent  was  to  find  a  partial  inversion  or 
partial  row  elimination  scheme  which  would  allow  for  the  solu¬ 
tion  of  the  boundary  probabilities  without  having  to  use  all 
the  information  in  the  matrix.  The  conclusion  is  that  the 
boundary  conditions  form  an  inseparable  portion  of  the  entire 
system.  Perhaps  further  research  may  uncover  an  efficient 
scheme  for  solution. 
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D.  AN  ATTEMPT  TO  SOLVE  FOR  THE  MARGINAL  PROBABILITIES 
USING  A  SEMI-MARKOV  MODEL 


> 

Ml 

M2 

Let  N^  =  #  of  customers  in  queue  #1 , 

Consider  just  the  N^  queue  (including  customer  in 
service) , 

N^  is  dependent  upon  the  arrival  stream,  the  service 
distribution,  and  the  distribution  of  N^.  Specifically, 
when  N2  =  N,  arrivals  cannot  enter  either  queue. 

The  N^  queue  slIso  has  capacity  =  N. 

Consider  modeling  the  N^  queue  as  a  semi-Markov  process. 

The  state  of  the  system  will  be  the  number  of  customers  in 
the  system.  N^=0,l,2 . N. 

Let  P.  .  be  the  transition  probability  from  i  — >  j . 

^  J 

Let  i  =  expected  time  system  is  in  state  i. 

Let  nr i  =  limiting  transition  probability  of  transitioning 
into  state  i. 

Let  P^  =  steady  state  probability  of  being  in  state  i. 

Transitions  will  only  be  allowed  from 

i— ^  i+1  and  i — »  i-1. 

Using  a  semi-Markov  model,  one  assumes  that  when  the  system 
Is  in  state  i  it  will  remain  in  state  i  for  a  random  amount 
of  time  (mean  =  i)  .  Then  the  system  will  either  transition 
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to  state  i+1  or  i-1.  When  1=0  it  will  transition  to  i=l  with 
probability  1.  When  i=N  transition  i=N-l  with  probability 

1. 

In  essence,  the  semi -Markov  model  has  transitions  which 
are  in  accordance  with  a  Random  Walk  Markov  Chain  (with  pure 
reflecting  barriers) . 

Thus  the  transitions  are  described  as 


0 

1  2 

3 

4 

1 — 1 

1 

N 

0 

0 

1  0 

0 

0 

0 

0 

1 

h 

0 

0 

0 

0 

0 

2 

0 

Pj  0 

^2 

0 

...  0 

0 

3 

0 

0  P3 

0 

^3 

...  0 

0 

4 

0 

• 

• 

• 

0 

N-1 

P  ,  0 

n-1 

^n-1 

N 

0 

0  0 

0 

0 

1 

0 

where 

Pi  +  li  = 

=  1 

The  underlying 

theory 

of  this  model 

allows 

the  following 

statements 

• 

If 

i  i=0,l,2. 

.  .  N 

are  the  solutions  to 

the  following 

system. 

N 

T  ■'Ti  =  1 

i=0 


N 

j=0 


31 
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Then  the 
Pi  =  ■ 


P.  ’s 


j=0 


can  he  found 
i  


ty 


The  basic  requirement  for  a  semi-Markov  model  is  that  the 
distribution  of  time  spent  in  a  state  be  independent  from  state 
to  state.  It  appears  as  though  the  exponential  character  of 
the  service  and  arrival  distributions  will  insure  this. 

It  should  be  noted  that  all  T^i’s  and  ’  s  must  be 
conditioned  on  the  state  of  the  N2  queue.  A  Bernoulli  type 
conditioning  will  be  used.  The  states  of  will  be  partitioned 
into  two  classes  N  P(N2  N)  =  1-K 

N  P(N2  =  N)  =  K 

1 .  Transition  Probabilities 
Consider  i  =  1,2,3»  •••  N-1 
Case  I  ^2^  ^ 

Given  the  system  is  in  state  i 

p  =  -A_ 
i ,  i+ 1  X  + 

p _ 

^i,i-l 

Case  II  =  N 

When  ^2  is  capacitated,  the  transition  from  i  to 
i+1  can  only  take  place  if  N2  serves  one  customer  and  then 
an  arrival  occurs  before  can  serve  one  customer. 


arrival  occurs  before  service 


service  occurs  before  arrival. 
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i.e.  S^'^EXP(-U^)  S2^EXP(M2)  A '^EXP(  \  ) 

* 

where  S^  represents  memoryless  service 


P.  . 


U 


\ 


\m 


i,i+l  M^+ 


(\  +A  ^)(  M^) 


Likewise  transition  i— >  i-1  takes  place  only  when  ^2"^^ 


=  P(S<rS2+A)  =  PCS^^S^)  +  P(S^>S2,  S^<A) 

JU  ^  ^  JU^{\  + 


P.  . 


^1  +  '^2 


i,i-l  '^l'^ 


'1  (  X  + '^^)  ( +  >u^) 


Thus  finding  the  total  transition  probabilities  as  a  function 
of  K,  one  has 

\  ,  .  X>Uo  K 

^i,i+l  ~X+i4^  (1-K)  +  (  \ 


h.i-i  =x-ni^  (i-K)  ^  (■xx-i:.^)ri7^ 


K 


2) 


Simplifying 


1 , 1-1 


p. 


^  ^{\k+  ^2) 

(  \+  (  'U^+ 

\  iM^+  1^2  -  '^l 


^  (  \+M^)  (  ^2) 


=q 


For  i  =  0 


^^0,1=  1 


i  =  N 


^N.N-1  ■  ^ 


/ 
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2.  Limiting  Transition  Probabilities 

Using  the  relationship  IT .  =  ^  P..  T  . 

1  j=0  ji  j 


the  following  system  of  equations  develop. 


r  = 

N. 

^  N-l 

nr 

=  '^  +0“'^ 

N-l 

N  ^  N-2 

=  qT  +  p  ‘tr 

M 

N-2 

^  N-3  _  N 

N-2 


T 

TT 

IT 


3 

2 

1 

0 


=  qT^  +  P 

=  q  Tr^  +  p 

s'Tr  +  p  Tr 
0  .  2 


Setting  the  equations  as  in  a  matriXj  one  row  operation 
through  the  matrix  produces  the  following  solutions 

Q  i=l,2,3,  ...N-1 

N  '^P  /  "  0 

Using  the  normality  relationship  ~  ^ 


IT  =  _ 2P^,1  . 

°  2(P-q(Jf-^7 
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3*  Expected  Time  in  Each  State. 


The  ^i's  are  calculated  by  conditioning  on 


Case  I  N2<  N 


i  =  0  ^  i  = 


1 


X 


i  =  1,2 


9  ••  • 


^  i 


Case  II  =  N 

i  =  0  the  queue  must  wait  for  1  service  and 
then  1  arrival 


i  =  1,2 . N-1  the  queue  waits  for  the  minimum 

of  1  service  or  1  N2  service  and  1  arrival. 

Let  X  =  min 


By  conditioning  etc. 


(  \  -  +  ^2) 


i=N 


c<  ^  _1 


still. 


i  li 


1 
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Hence  the  total  expected  waiting  times  are 


/\k+  >U 


.U 


X  K+  H  -tl  2 


i  =  0 


rx+ u^) 


i  =  1,2,.. .N-1 


1 

>u-, 


i  =  N 


Knowing  and  oC^'s  one  can  calculate  closed  form  expres¬ 
sions  for  in  terms  of  \  , -U^,  and  K. 

Now  from  original  two  dimensional  steady  state 
transition  diagram  one  can  show  the  following  relationship  to 
be  true.  This  equation  was  derived  in  the  previous  section. 


(  X  +  M  ^ 


^1  =^^1  ^2  ^ 


0 


For  the  case,  K=3.  the  above  equation  could  not  be 
satisfied  using  the  expressions  for  found  by  the 

proposed  model.  It  appears  as  though  the  fixed  coupling 
coefficient  K  does  not  completely  characterize  the  interaction 
between  the  two  queues.  At  least  having  the  above  marginal 
equation  available  allowed  for  a  quick  check  against  the 
model . 


54 


V .  CONCLUSION 


There  is  a  great  deal  more  research  which  can  be  performed 
in  order  to  better  imderstand  the  effects  of  arrival  dependence 
upon  service  systems. 

As  one  can  no  doubt  conclude,  finding  analytical  solutions 
to  queueing  problems  of  this  type  is  extremely  difficult  and 
may  in  fact  require  a  completely  novel  approach. 

The  infinite  server  queue  did  not  pose  much  difficulty 
because  the  waiting  time  of  any  customer  was  strictly  a  fionction 
of  his  own  service  time.  However,  the  systems  with  finite 
servers  and  finite  capacity  still  remain  unsolvable  from  an 
analytical  viewpoint. 

One  can  think  of  many  instances  where  finite  server  and 
finite  capacity  queues  could  be  used  to  model  repair  and 
maintenance  systems.  Consider  a  military  maintenance  company 
which  services  vehicles  or  generators.  Often,  the  arrival 
of  a  certain  item  of  equipment  initiates  several  activities, 
many  of  which  are  independent.  It  was  the  aim  of  this 
research  to  model  such  activities.  Unfortunately,  only 
partial  answers  have  been  obtained.  Hopefully,  more  study 
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will  be  concentrated  in  this  area  of  stochastic  modeling  in 
order  to  advance  the  current  knowledge  of  these  and  other 
related  systems. 
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